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PREFACE 


This  report  was  prepared  by  Pacif ic-Sierra  Research  Corporation 
(PSR) ,  a  subsidiary  of  Eaton  Corporation,  for  the  Defense  Nuclear 
Agency  (DNA) .  The  work  was  done  under  contract  DNA  001-85-C-0089,  and 
supervised  by  Dr.  Michael  J.  Frankel. 

In  previous  PSR  studies,  the  generation  of  fire  winds  by  large 
area  fires  has  been  considered.  In  this  report,  the  possible 
mechanisms  for  generating  central  column  swirl  are  examined,  and 
induced  tangential  velocities  are  calculated  for  a  10-km  radius  large 
area  fire. 
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SECTION  1 
INTRODUCTION 


Large  area  fires,  sometimes  called  firestorms,  differ  in  many 
ways  from  small  concentrated  flames,  forest  fires  or  even  isolated 
urban  fires.  The  distribution  of  buoyancy  over  a  large  area  can 
produce  high  winds  that  flow  inward  from  all  directions  once  the 
firestorm  has  been  established.  The  fire  winds  are  the  source  of 
oxygen  for  the  fire,  and  also  tend  to  limit  its  spread  by  forcing  the 
flames  inward.  Since  the  burning  zone  does  not  move,  another 
phenomenon  typical  of  firestorms — complete  burning  of  all 
combustibles — can  result.  From  dimensional  arguments,  the  average 
inflow  velocity  at  the  fire  parameter  can  be  shown  to  grow  with  the 
cube  root  of  the  equivalent  fire  radius.  Thus,  only  fires  of  area 
above  3ome  minimum  will  necessarily  be  stationary  and  develop  into 
firestorms.  That  area  (or  typical  radius)  is  obtained  when  the  inflow 
velocity  is  higher  than  the  greatest  wind  a  fire  can  spread  against. 
That  value,  although  not  well  established,  is  evidently  a  function  of 
the  density  of  available  combustibles,  the  distribution  of  ignitions, 
and  the  shape  of  the  overall  fire  front.  One  value  quoted  in  the 
literature  [Glasstone  and  Dolan,  1977]  is  8  mph  or  3.5  m/s,  but  that 
value  is  probably  low,  as  indicated  by  the  February  19^5  Dresden  fire 
[Irving,  1963].  A  fire  diameter  of  0 ( 1 0 3)  m  is  probably  sufficient  to 
be  stable  in  such  ambient  winds.  Area  fires  of  that  size  are  unusual 
which  may  explain  the  rarity  of  firestorms.  Nevertheless,  natural 
disasters  such  as  volcanic  eruptions  or  earthquakes  could  conceivably 
Ignite  large  areas  simultaneously  so  that  large  inflow  velocities  are 
obtained . 

Large  area  fires  have  been  produced  artificially  by  patterned 
simultaneous  ignitions  in  firebombing  during  World  War  II,  by  planned 
ignitions  for  forest  clearing  [Sma^l,  1985],  and  by  nuclear  explo¬ 
sions. 
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It  has  been  suggested  [Emmons,  1965;  Long,  1967]  that  the  high 
winds  and  exceptionally  tall  plumes  related  to  firestorms  require  the 
existence  of  swirling  motions.  Sketchy  (mostly  anecdotal  eyewitness) 
evidence  from  the  19^3  Hamburg,  Germany  fire  has  been  since  quoted 
[Carrier  et  al.,  1985]  to  support  that  notion.  If  true,  the  existence 
of  swirl  has  critical  bearing  on  the  plume  generated  by  a  large  area 
fire  and  as  a  consequence,  on  the  injection  and  spread  of  combustion 
products  through  the  atmosphere  and  on  the  possibility  of  a  nuclear 
winter.  Massive  swirling  hinders  entrainment  and  mixing,  narrows  the 
plume  (by  centripetal  acceleration  of  the  lighter  plume  elements),  and 
drives  it  higher — possibly  into  the  stratosphere  where  particles  are 
likely  to  have  much  longer  residence  times. 


SECTION  2 

POSSIBLE  CAUSES  OF  SWIRL 


In  this  section,  we  discuss  the  possible  physical  mechanisms  that 
can  cause  a  large  area  fire  plume  to  swirl.  The  fires  considered  have 
burning  areas  greater  than  1  km2. 

Basically,  the  buoyancy  generated  by  an  area  fire  induces  a 
radial  inflow  of  ambient  air  at  ground  level  and  an  upward  convection 
of  combustion  products  [Larson  and  Small,  1982;  Small,  Remetch,  and 
Brode,  1984;  Small  and  Larson,  1986],  Thus,  swirl  is  not  an  intrin¬ 
sic  characteristic  [Morton,  1970]  of  the  flow  field  and  must  be  caused 
by  some  outside  agent,  for  example,  the  earth's  rotation.  The  area 
affected  by  the  fire  can  have  a  radius  of  0(10^)  m  so  that  Coriolis- 
type  effects  are  not  obviously  negligible. 

The  Coriolis  force  (acceleration)  can  be  written  as 

a  »  2u>  x  u  ,  ( 1 ) 

c 

where  w  is  the  earth's  angular  velocity  (actually  its  component  per¬ 
pendicular  to  the  surface,  w  -  7.27.1 0"5  sin#  rad/s).  The  inflow 
ve.ocity  u,  is  a  function  of  the  distance  from  the  fire  center.  The 
tangential  velocity  at  the  fire  rim  R0  at  time  T  due  to  the  Coriolis 
effect  is 


2uu  dt  »  - 


where  s  is  the  initial  distance  of  the  particle  that  reaches  the  rim 
at  time  T.  A  mature  flow  field  is  assumed  at  T  «  0.  In  reality,  the 
flow  develops  over  some  period  of  time.  The  minus  sign  indicates 
inward  velocity  u.  That  velocity  will  increase  with  time,  a  phenomenon 


R 

f 


2ii)U  dr/u 


2u)(s  -  R0) 


(2) 


I, 
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usually  called  "spin  up."  However,  the  fire  only  burns  for  a  limited 
time,  which  can  be  written  as 


t 


max 


me 

J 


(3) 


where  m  is  the  combustible  loading  mass  per  unit  area,  C  the  caloric 
equivalent  of  the  combustibles  per  unit  mass,  and  J  the  rate  of  con¬ 
version  of  chemical  energy  (burning)  per  unit  area.  The  latter  is 
sometimes  presented  per  unit  volume,  in  which  case  the  flame  height  Hf 
needs  to  be  considered  (i.e.,  (J)area  “  (^volume  *  Hf)*  The  elapsed 
time  t,  can  be  written 


where  the  inflow  velocity  u  was  assumed  similar  to  that  generated  by  a 
sink  flow  [i.e.,  inversely  dependent  on  distance  (u  -  1/r)].  This 
implies  that  the  inflow  into  a  large  fire  system  is  mostly  contained 
in  a  shallow  surface  layer.  Such  flows  have  been  calculated  [Small, 
Remetch,  and  Brode,  1984;  Penner,  1985].  From  Eqs.  (2)  and  (4),  and 
recalling  that  t  £  traax , 


Using  typical  combustible  loadings  (see  Table  1 )  and  calculations 
of  large  area  fire  flow  fields  as  input  conditions,  we  obtain 
(v/U0)max  of  approximately  0.4  to  0.5  (for  a  circular  area  fire, 
toward  the  end  of  the  burn).  That  is,  the  total  velocity  is  never  more 
than  13  percent  above  the  radial  influx  alone.  Data  for  the  Hamburg 


fire,  which  has  been  quoted  as  possibly  swirling,  has  been  included  as 
case  III  in  Table  1.  It  is  clear  that  some  swirl  can  be  generated  by 
the  earth's  rotation,  but  it  is  much  too  small  to  contribute  sig¬ 
nificantly  to  the  velocity  field.  A  typical  hurricane  calculation 
(case  IV)  shows  that  long  spin  up  times  are  required  to  cause  very 
large  increases  in  the  wind  speeds.  The  results  of  those  simple  cal¬ 
culations  are  consistent  with  observations.  Thus,  it  3eems  reasonable 
that  the  present  model  properly  estimates  the  potential  for  swirl  from 
Coriolis  accelerations. 

A  different  possible  source  of  swirling  motions  lies  in  non¬ 
uniformities  in  the  burning  area,  or  its  surroundings.  Those  in¬ 
homogeneities  can  be  in  the  type  and  concentration  of  combustibles,  in 
the  topography  of  the  area,  or  in  the  time  interval  of  the  ignitions. 
Each  of  those  effects  can  cause  nonradial  (tangential)  components  in 
the  flow  field.  Conservation  of  angular  momentum  will  cause  the  tan¬ 
gential  motions  to  be  intensified  as  the  flow  elements  involved  move 
radially  closer  to  the  fire  epicenter. 

To  describe  that  effect  in  a  more  quantitative  manner,  consider 
an  example  of  a  square  city  block  of  side  length  a,  which  is  burning 
uniformly,  but  the  access  routes  for  incoming  air  (the  surrounding 
streets)  are  blocked  (e.g.,  by  rubble).  All  the  air  feeding  the  fire 
comes  in  from  two  channels  (streets)  of  width  b,  such  that  the  vor- 
ticity  is  increased  rather  than  reduced.  That  is  a  worst-case 
scenario,  as  any  air  coming  in  from  other  streets  will  reduce  the 
asymmetry . 

The  circulation  induced  by  the  nonradially  centered  inflow  is 


T  *  2u1a 


(6) 


where  u1  is  an  appropriately  averaged  speed  in  the  street.  This  can  be 
written  as 


5 


Table  1 .  Calculation  of  swirling  component  of  velocity  at  ground  level  in  vicinity  of 
large  area  fire. 
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Time  to  total  burnout  (actual  full  strength  fire  duration  for  case  III). 
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where  M  is  the  mass  of  air  drawn  in  per  unit  time,  p  is  its  density, 
and  h  the  inflow  region  height.  The  air  flux  can  be  related  to  the 
fire  by 


where  B  is  the  fuel  loading,  S  the  stoichiometric  ratio,  L  the  air 
excess  ratio,  and  t  the  burning  time,  weighted  so  that  a  uniform 
burning  rate  is  assumed.  From  Eqs.  (6)  through  (8) 


LSBa  _  Ma 
phtb  “  phb  ’ 


i.e.,  the  circulation  produced  is  proportional  to  the  burning  rate  and 
burning  unit  size  and  inversely  proportional  to  the  available  air 
entry  channel  size. 

Each  of  those  variations  can  happen  on  different  scales.  For 
example,  individual  buildings  on  a  block  may  burn  at  different  rates, 
or  produce  different  levels  of  heat  release  per  unit  time.  This  is  an 
0(10)-m  scale.  Streets  will  usually  not  burn,  so  that  blocks  are 
separated.  Next,  there  could  be  a  gradation  of  buildings  from  residen¬ 
tial  to  commercial  and  industrial  zones,  i.e.  a  scale  of  hundreds  of 
meters  to  kilometers.  Any  variations  that  have  length  scales  much 
smaller  than  the  equivalent  radius  of  the  fire  are  probably  not  sig¬ 
nificant,  as  their  effects  on  rotation  will  cancel  out  on  the  3cale  of 
the  entire  fire.  Thus,  if  a  typical  block-edge'  causes  swirling  in  the 
clockwise  direction,  there  is  a  high  probability  that  opposite  rota¬ 
tion  of  similar  strength  will  occur  in  another  street.  The  total  over 
many  blocks  should  roughly  cancel  out  with  the  residual  circulation 
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being  inversely  proportional  to  the  number  of  units  of  combustion  N 
(burning  blocks).  That  can  be  written  as 


r.  .  ,  oc 
total 


(10) 


The  sum  tends  to  zero  as  N  becomes  large. 

Alternatively,  variations  on  scales  much  larger  (hundreds  of 
kilometers)  than  the  fire  scale  are  probably  negligible  too  as  the 
tangential  components  of  flow  they  produce  are  very  small. 

Thus,  we  are  left  with  nonuniformities  of  the  3cale  of  the  fire. 
Those  might  be  topographical  (a  hill  or  mountain,  or  a  lake)  or  a 
change  in  neighborhood  character  as  described  previously.  Before 
discussing  them  quantitatively,  further  details  of  the  relation  be¬ 
tween  those  nonuniformities  and  swirling  motions  are  needed. 

A  nonradial  inflow  (streamline)  pattern  can  be  expected  if  the 
burning  area  is  not  axisymmetric ,  since  the  local  (near-field)  inflow 
is  normal  to  the  fire  front.  Thus,  there  will  be  spreading  of  stream¬ 
lines  near  concave  parts  of  the  perimeter,  and  concentrations  next  to 
convex  sections.  Perimeter  shape  effects  alone,  however,  will  not 
produce  large  scale  net  vorticity,  as  the  tangential  components  cancel 
out  for  a  uniform  strength  fire  in  a  closed,  simply  connected  region 
(i.e.,  r  -  f  V  •  ds  -  0).  The  cancellation  may  not  be  perfect, 
however,  if  the  burning  is  nonuniform.  For  example,  an  elliptical  or 
lunar  shaped  area  with  greater  burning  intensity  at  one  end  can 
produce  a  net  tangential  component  at  any  given  circle  surrounding  the 
fire  zone.  Also,  if  a  nonuniform  burning  area  is  located  (for  example) 
on  the  side  of  a  hilly  ridge,  the  horizontal  inflow  asymmetry  is 
enhanced,  leading  to  possible  net  vorticity. 

Detailed  predictions  of  the  possible  swirling  motions,  produced 
by  nonuniformities  in  the  fire  and  f.uel  bed,  depend  on  a  great  variety 
of  possibilities.  Some  estimates  of  the  upper  limits  of  rotation 
resulting  from  this  class  of  vorticity  generators  are  presented. 
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We  now  estimate  the  vorticity  produced  by  a  10-km  radius  fire. 

The  example  city  is  bordered  on  one  side  by  a  steep  mountain  range. 

The  fire  burns  an  asymmetrically  loaded  area,  with  average  combustible 
density  on  the  right  half  of  the  fire  zone  twice  that  of  the  left. 
Similarly,  the  burning  rate  is  twice  as  large  in  the  right  half.  The 
heat  release  QA  in  the  left  fire  zone  is  assumed  1  kW/m^;  roughly  that 
from  an  area  fire  in  an  average  density  city.  A  uniform  heat  release 
of  that  magnitude  can  generate  inflow  speeds  of  about  13  m/s  in  the 
"mature"  stage  of  the  fire  [Small,  Remetch,  and  Brode,  1984].  We  now 
assume  that  the  mountains  reduce  the  inflow  from  the  sector  they 
surround  by  a  factor  of  2  (over  one-half  of  the  circumference).  Each 
of  these  inhomogeneities  (topography  and  fire  intensity)  separately 
would  not  cause  swirling,  but  the  combination  conceivably  could.  The 
offset  distance  A  of  the  weighted  center  (from  the  geometrical  center) 
of  the  fire  is  related  to  the  heat  release  by 
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In  this  example,  Qg  «  2  kW/m^ ,  QA  -  1  kW/m3 ,  and  A  »  0.141  R  -  1.4  km. 
Assuming  a  linear  relationship  between  heat  release  rate  and  induced 
velocity,  then  Vg  -  2  VA.  This  is  not  strictly  true  [Small  and  Larson, 
1986],  but  is  a  sufficient  approximation  for  this  calculation. 

The  circulation  Induced  by  the  asymmetries  in  burning  and  topog¬ 


raphy  is 
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The  constant  c  relates  the  heat  release  to  the  induced  fire  winds. 
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The  circulation  is  then  r  •  3-8  x  to1*  n^/s,  and  the  induced  tangential 
velocity  (at  an  altitude  above  the  hills)  about  0.60  m/s.  This  assumes 
no  necking  of  the  plume.  If  the  plume  cross-section  area  is  reduced  by 
3/4  (half  the  fire  radius),  then  the  Induced  tangential  velocity  is 
only  1.2  m/s.  For  either  a  uniform  heat  release  or  uniform  topography, 
no  swirl  component  is  introduced.  For  the  combination  of  asymmetries 
and  even  with  an  extreme  necking  of  the  plume  (worst-case  analysis), 
the  swirl  generated  is  essentially  negligible. 

Perhaps  the  mechanism  most  likely  to  induce  swirling  [Carrier  et 
al.,  1985]  is  vorticity  already  existing  in  the  atmosphere — brought 
in,  and  concentrated  by  the  ground  level  inflow  of  the  fire.  We  can 
reasonably  exclude  as  a  source  extreme  weather  systems  such  as  3torms 
or  hurricanes  (1)  because  of  their  relative  rarity,  and  (2)  because 
they  will  probably  help  subdue  the  fires,  and  in  any  case,  cause  much 
larger  changes  in  any  plume  before  merging  with  it  to  produce  a  swirl¬ 
ing  plume.  The  effect  of  small  eddies  [R  -  0(100)  m  or  less]  can  also 
be  neglected;  they  occur  frequently,  but  their  effect  cancels  statis¬ 
tically  or  quickly  decays  with  time  (at  a  rate  Inversely  proportional 
to  the  eddy  size  and  viscosity). 

The  greatest  influence  may  result  from  horizontal  wind  gradients 
or  eddies  of  wavelength  comparable  to  the  fire  diameter.  Such 
gradients  and  eddies  can  be  produced  by  winds  passing  suitable  topog¬ 
raphy  upstream  of  the  fire  [Hunt  et  al.,  1978].  The  effect  is  equiv¬ 
alent  to  having  a  tangential  (swirl)  velocity  of  one-half  the  field 
velocity  difference. 

A  reasonable  upper  limit  for  velocity  differences  from  such 
gradients  is  about  10  m/s,  which  means  that  a  tangential  component  of 
5  m/s  1 s  expected.  Such  values  would  be  characteristic  of  windy  condi¬ 
tions.  Higher  values  could  result  under  storm  conditions.  Neverthe¬ 
less,  such  shears  are  not  large  enough  to  cause  rapid  swirling.  The 
vorticity  generated  is  of  similar  order  as  those  previously  discussed. 

In  view  of  the  above,  It  seems  justified  to  reexamine  popular 
conceptions  about  the  Hamburg  fire.  The  positioning  of  the 
eyewitnesses — especially  those  in  streets  relatively  close  to  the  fire 
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perimeter,  in  part  explains  the  reported  impressions  of  "very  high 
winds  and  swirling."  The  inflow  winds  are  in  general  channeled  down 
streets  and  probably  more  so  in  old,  tightly  packed  cities  like  Ham¬ 
burg.  That  causes  an  increase  in  speed  due  to  two  reasons:  (1)  the 
decrease  in  available  area,  which  may  be  a  factor  of  2  to  3;  and  (2) 
the  Borda  effect  [Milne-Thomson,  1968]  which  causes  flow  in  con¬ 
stricted  areas  to  neck  down  further  and  utilize  only  about  one-half  of 
the  available  area.  Thus,  for  city  streets  or  natural  channels,  in¬ 
creases  by  a  factor  of  4  to  6  can  be  expected  over  the  global  or 
average  inflow  velocities  (see  Table  1).  Accordingly,  instead  of  7  m/s 
(case  III  of  Table  1)  it  may  not  be  surprising  to  observe  peak  radial 
speeds  of  35  m/s  ( ~  80  mph) . 

Also,  the  streets  were  not  necessarily  perpendicular  to  the  flame 
perimeter  so  that  flow  down  those  urban  canyons  could  be  interpreted 
as  having  a  tangential  (swirling)  component.  Local  fire  whirls  could 
be  shed  at  various  locations  along  the  perimeter.  However,  the  total 
vorticity  for  the  whole  fire  system  is  still  negligible.  Thus  the 
impressions  of  ground-level  observers  do  not  necessarily  infer  large- 
scale  organized  plume  swirling  motions. 

Finally,  swirl  is  not  a  prerequisite  to  the  generation  of  large 
velocity  fire  winds,  Larson  and  Small  [1982a, b]  and  Small  and  Larson 
[1986]  have  demonstrated  analytically,  and  Smith,  Morton,  and 
Leslie  [1975],  Small,  Remetch,  and  Brode  [198U],  and  Penner  [1985], 
have  shown  in  hydrocode  simulations  that  pressure  gradients  generated 
by  the  f ire-produced  buoyancy  distribution  account  for  the  high- 
velocity  fire  winds. 
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SECTION  3 
CONCLUSIONS 


The  analysis  presented  in  the  previous  section  shows  that  all  the 
factors  considered  do  not  cause  appreciable  swirling  by  themselves, 
even  though  the  most  extreme  conditions  were  examined.  Thus,  for  the 
question,  Can  the  plume  resulting  from  a  large  area  fire  swirl?  the 
answer  is,  probably  not.  For  rather  unusual  combinations  of  condi¬ 
tions,  such  as  nonuniform  fires  burning  for  a  long  period  (very  high 
combustible  loading  not  easily  accessible  to  air) ,  in  an  area  par¬ 
tially  bounded  by  hills,  swirl  velocities  large  enough  to  be  notice¬ 
able  (although  not  large  enough  to  significantly  affect  the  plume  rise 
in  the  atmosphere)  may  occur — especially  at  the  later  stages  of  the 
fire.  That  explains  why  most  observed  large  fire  plumes  [including 
experimental  fires  such  as  Meteotron  (Church,  Snow,  and  Dessens, 

1980),  Flambeau  (Countryman,  1969),  and  Chapleau  (Small,  1985)]  did 
not  swirl. 
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